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1. Introduction 

The Woods-Saxon potential and its various modifications have been received for much 
interest describing metallic clusters in a successful way. It is used in the central part 
of the interaction neutron with one heavy-ion nucleus and also for the optical potential 
model PJ. Differential cross-section in the ^^O C elastic scattering are analyzed 
in some energies by using this potential |2]. The quantum behavior of the relativistic 
particle in the presence of a central Woods-Saxon potential has settled the possible 
existence of bound state spectra. The relativistic Dirac-oscillator and Dirac-exponential- 
type potential problems have already been established by adding an off-diagonal linear 
radial term to the Dirac operator for a long time ago [Sllll- Recently, the relativistic 
bound states spectrum and its eigenfunctions for the triaxial and axially deformed 
harmonic oscillators have been derived as well jS]. In addition, a mixture of the Dirac 
oscillator (tensor potential) with vector and scalar harmonic oscillator potentials has 
been solved analytically for the general case [E]. 

Furthermore, only a few articles for the relativistic problems have been written on 
the Dirac equation with the exponential-type potential. The Dirac equation has been 
solved by making use of two-component spinors for the exponential type potentials such 
as Woods-Saxon and Hulthen potentials for a special case. Kennedy has studied the 
generalized approach to the Woods-Saxon potential and obtained the scattering and 
bound-state solutions of the one- dimensional Dirac equation. However, more realistic 
cases have not been discussed [Zj. A. D. Alhaidari has just introduced a new formalism 
to the definition of the radial Dirac equation and solved for a class of shape-invariant 
potentials |H1 IHl ■ The main point in the formalism is that two coupled first order 
differential equations resulting from the radial Dirac equation generate Schrodinger-like 
equations for the two spinor components. Following the procedure given in Ref. [S], 
we present a new systematical approach to solve the Dirac- Woods-Saxon problem by 
means of the Nikiforov-Uvarov (NU) method JT]. The non- relativistic limit reproduces 
the well-known non-relativistic energy spectrum and results in the Schrodinger equation 
for a special form of the Woods-Saxon potential [12]. 

The article is structured as follows: In Section 2, we briefly introduce an overview 
of the technical details of the formalism improved by Alhaidari. After that, the basic 
concepts of the Nikiforov-Uvarov method are given in the same section to solve the 
Dirac- Woods-Saxon problem. Section 3 is devoted to the solution of the problem to 
obtain the energy eigenvalues and eigenfunctions by applying the NU method. The 
paper is concluded with a short summary in Section 4. 

2. Formalism and Method 

We shall first introduce the Alhaidari's formalism proposed to solve the Dirac equation 
for spherically symmetric potential interactions. Later, the fundamental mathematical 
relations of the NU method will be summarized to obtain the solution of the Schrodinger- 
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like equations easily and systematically. 



2.1. Overview of the Formalism 

Alhaidari's formalism is based on the consideration by writing the relativistic 
Hamiltonian for a Dirac spinor coupled to a four-component electromagnetic potential 
{Aq,A). After using Gauge invariance and spherically symmetric of the electrostatic 
potential, free Dirac equation transforms to the matrix representation of the Dirac 
Hamiltonian (see Ref.[8j for more detail). For convenience, atomic units are selected as 
m = e = h = 1 and the speed of light c is matched with a~^. Thus, the Hamiltonian 
for a Dirac spinor in four-component electromagnetic potential {Aq, A) can be written 
as follows: 

/ 1 + aAo —iaa.V + iao.A \ 

(1) 

—iaa.'V — iaa.A —1 + aA^ y 

where a is the fine structure parameter and a are the three 2x2 Pauli spin matrices. 
Taking the spherically symmetric case and writing {Aq,A) as {aV{r).,fW{r)), the two- 
component Dirac equation is obtained as 

/ l + a^V{r)- Er 



H 



\ 



H 



a 



a 







(2) 



y + W{r) + j-^\ -l+a^V{r)-Eji J 

where /(r) and g{r) are real radial square integrable functions, Er is the relativistic 
energy and k, is the spin-orbit coupling parameter defined as/t = ±(j + l/2) = ±1, ±2, ... 
for / = j ± 1/2. However, the current problem is analytically solvable only for £ = (s- 
states). In addition, V{r) and W{r) are the even and odd components of the relativistic 
potential, respectively. For a given value of the spin-orbit coupling parameter n, 
Schrodinger-like requirement relates the two potential functions as: W{r) = ^V{r) — ^, 
where is a real parameter and V{r) is not depend on the k parameter. In order 
to obtain the Schrodinger-like equation in the formalism proposed by A. D. Alhaidari, 
it is used a global unitary transformation which eliminates the first derivative. Thus 
U{r)) = exp{^arja2) is applied in Eq.Q. Here, ?7 is a real constant and a2 is the 2x2 
Pauli matrix which defines the two radial spinor components in terms of the other. 



a 



C±Er 



0±(r), 



with C 



cos^arf) = ^1 — (aC,)'^ > 0, 




(3) 



(4) 



Here 0^(r) is the upper or lower spinor components respectively. It is emphasized that 
Eq.Q with the top and bottom signs are not valid for negative and positive energy 
solutions respectively. The top and bottom signs in front of Eji in Eq.Q are not 
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allowed to take the values —C and +C. Because these two values are elements of the 
negative and positive energy spectra respectively. Substituting these into the radial 
Dirac equation (Eq.Q), we get the Schrodinger-like second-order differential equation 
in terms of the lower and upper spinor components as 



C'2 . X. CdV El 



dr- ' e- ' — ' idr a- ^^^'^ = ^^ 

where the "+" sign belongs to the upper spinor component, while the other sign 
corresponds to the lower one. 

2.2. Basic Concepts of the Method 

Solutions of the Schrodinger-like second order differential equations play an essential 
role in studying many important problems of theoretical physics. In this point, the 
NU method can be used to solve these types equations with an appropriate coordinate 
transformation s = s{r) [TT] : 

n-s) + ^^fis) + ^Ms) = (6) 

where a{s) and a (s) are polynomials with at most second-degree, and r (s) is a first- 
degree polynomial. It is of fundamental importance in the study of particular special 
orthogonal polynomials ^Hl- These polynomials try to reduced Eq.® to a simple form 
by taking ip^s) = (j){s)y{s) and choosing an appropriate Consequently, Eq.(jni) can 

be reduced to an equation of hypergeometric type 

a{s)y'\s)+T{s)y\s) + Xy{s) = 0, (7) 

where r(s) = r (s) -|-27r(s) (its derivative must be negative) and A is a constant, which 
is given in the form 

A = A„ = -nr' - ^^^^^a", (n = 0, 1, 2, ...). (8) 

Here, A or A^ are obtained from a particular solution of the form y{s) = yn{s) which is 
a polynomial of degree n. yn{s) is the hypergeometric type function whose polynomial 
solutions are given by Rodrigues relation 

where i?„ is the normalization constant and the weight function p{s) must be satisfied 
the condition 

[a{s)p{s)]' = T{s)p{s). (10) 

To determine the weight function given in Eq. ljlOj) . we must immediately obtain the 
polynomial 7r(s) from: 



7T = ± 



^ 2 

a'— T ^ 



a+ka. fill 
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In principle, the expression under the square root sign in Eq. ljlip can be arranged 
as the square of a polynomiaL This is possible only if its discriminant is zero. In 
this case, it is obtained an equation for k. After solving this equation, the obtained 
values of k are included to the NU method and here there is a relationship with A of 
k so that k = X — 7r'(s). After this point, an appropriate can be invented from 
0(s)V0(s) = 7r{s)/a{s). 



3. Special form of the Woods-Saxon potential 



The interaction between nuclei is commonly described by using a potential that consists 
of the Coulomb and the nuclear potentials. It is usually taken in the form of the Woods- 
Saxon potential. Here, we take the following special form for the Woods-Saxon potential 
which is specified by "g" parameter 

qVo 



V{r) 



(12) 



where Vq is the potential depth, Rq is the width of the potential, b is thickness of the 
surface which is usually adjusted to the experimental values of ionization energies and q 
is a real positive parameter which is responsible for the specification of the Woods-Saxon 
potential. After substituting the potential into Eq.(jSl), we obtain an equation for the 
upper spinor component 




2qERVo qCVo 



El 



1 



q + 



qCVo —f^ 



2qERVo 



0+(r) = 0, 



El 



q + 



q + 



0+(r) = 0. (13) 



In order to apply the NU— method, we rewrite Eq.(|13|) by using a new variable of the 
form s = —e^^~b^\ 



s d ( s d 
h ds \ h ds 



q'^C^Vi 



1 — qs 



^ 2qsEjiVo ^ qCVo 



El 



1 



By introducing the following dimensionless parameters 



'El 



l-qs ■ ib {l-qsf 
qCVob 



P = 2qERVob' 



7 



0+(s) = 0.(14) 



(15) 



we reach the following hypergeometric type equation defined in Eq.( 



(s)^ l-qs 



+ 



■X 



{q^e + ql3- 7 )s - i2qe + f3 + -f)s + e 



ds"^ s{l — qs) ds s^(l — qs)'^ 
After comparing Eq.()16p with Eq.®, we obtain the corresponding polynomials: 

r(s) = l-gs, a{s) = s{l-qs), a (s) = {q'^e+qP-j^)s^-{2qe+P+j)s+e.{17) 



0.(16) 
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Substituting these polynomials into Eq. (fTT|) . we organize the polynomial n{s) as follows 

= -y ± ^ - + - - 4^^) + 4 (2^9 + /5 + 7 + A;) s - 4e, (18) 

with cr'(s) = 1 — 2gs. It is taken into consideration that the discriminant of the second 
order equation under the square root sign has to be zero. Hence, the expected roots are 
obtained as k± = {—[3 — 7 ± z(g + 27)1/5). In this case, substituting these values for 
each k into Eq.(fTH|). the possible solutions are obtained for 7r(s) 

(g + 27 - 2iq^) s + 2iy/^, 
for /c+ = -/5 - 7 + i(g + 27)1/5 

(19) 

(g + 27 + 2iq^) s - 2i^, 
for k_ = —[3 — 'y — i{q + 27)1/5 

From the four possible forms of the polynomial vr(s), we take a certain one which is 
the derivative of r(s) has a negative value. Therefore, the function r(s) satisfies the 
following equalities: 

r(s) = 1 + 2iy/6 - s (3g + 27 + 2iqy/e 



7r(s) 



1 

2 2 



r (s) = -(3g + 27 + 2zgv^) . (20) 
In the present case 

As) = - i [(g + 27 + 2zgv^ ) s - 2iq^e ] . (21) 

From k = X — vr'(s) and also Eq.(jHJ, we obtain respectively: 

A = -/3-27-g-2i(7 + g)v^ (22) 

A = A„ = n'^q + 2nq + 2^7 + 2niq^. (23) 

After having the comparison of Eq.(j221) and Eq.(j2Sl) and substituting the values of e and 
we can immediately obtain the K-independent relativistic energy eigenvalues Ejinq of 
the Dirac particle as follows 

^1, = -[2&(2^' + h\^Via^)\-\{\ + n)qV^h{J + 7)^2 ± + ^^2(3. + 7)'9'^o&' 

(24) 

_(T2 + 6VK'«')((1 + nfij + 7)2 - 462T2)]V2}^ 

where T = (1 + n)g + 7. To have a physical result, the expression under the square root 
must be positive, n is a positive integer defined in the interval of n^ax > > and is 
called the radial quantum number. 

By interesting with Eq. p3p . we can easily show that, in the non-relativistic limit 
a — i> 0, the relativistic energy is a limit of the non-relativistic energy, Eji fii 1 + a'^Ejyfi, 
where E^r is the non-relativistic energy. The wave equation is reduced to the following 
form, choosing g = 1: 

2Emr dy-^ir) = 0. 



7 

h — 



7 



1 + ey b . 



1 + e 



(^) 



(25) 
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To obtain a more suitable case, we can use the following form after taking 7 



Vo - 1/262 



1 + e^ 



2E 



NR 



0+(r) = 0, 



(26) 



which is in the form of the Schrodinger equation for the new type s-wave non-relativistic 
Woods-Saxon potential. The corresponding energy spectrum has already been given in 
Ref. P3] by using the hypergeometric functions, and also repeated by means of the NU 
method in Ref. ^H] as: 



E 



NR 



b{Vo - 1/262) n + 1 



+ 



(27) 



n + 1 ' 2b 

where the transformation Vq Vq — 1/26^ is applied for the convenience and the index 
n relates to the radial quantum number (n = 0,1,2,3,...). In order to obtain the 
relativistic energy spectrum directly, considering the relativistic Eq. ()13|) and the non- 
relativistic Eq.(j2Z|) for the case of 7 = —1, we can propose the relevant parameter 
map: 

6 — > 6, Ro ^ Ro, 

Vo - 1/262 _^ _ 1/262, 



E 



NR 



(28) 



{El - l)/2a\ 

Using the map between the parameters of the two equations, the resulting upper 
relativistic energy spectrum is found as follows: 



E 



-Voha^{{n + 1)2 - 1) + J462(n + 1)^ + {n + iya^[Ah*Vi - ((n -f 1) 



R 



2b{{n + 1)2 + 62^02^2) • (2^) 

It can be easily seen that the relativistic energy spectrum in Eq. ()24|) gives the same 
result as Eq. but it can be used only for the s-states. In addition, Eq. indicates 
that one deals with a family of the Woods-Saxon potential and the relativistic energy 
spectrum will be also used to describe the single-particle motion in nuclei. Therefore, 
Eq.(j21I) can also give the solution of the relativistic Dirac-Woods- Saxon problem with 
a general value of 7 for q = 1. 

Let us now find the corresponding wave functions. According to the NU- 
method, the polynomial solutions of the hypergeometric function y{s) depend on the 
determination of weight function p(s) satisfying the differential equation [o"(s)p(s)]' = 
t{s)p{s). Thus, p{s) is calculated as 

p{s) = (1 - qsY s''-^, (30) 

where u = 1 + Substituting into the Rodrigues relation given in Eq.Q, the wave 
functions are obtained in the following form 



Vngis) = v4„ (1 - qsY 



ds"' 



[I - qs 



(31) 



where An is the normalization constant. Taking q = 1, the polynomial solutions of 
yn{s) are expressed in terms of the Jacobi Polynomials, which is one of the orthogonal 
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polynomials. In this case, the weight function is (1 — s)'^s^*^ and Eq.(|31|) is reduced 
to ~ P^^'^'^^l - 2s) [ini. After substituting n{s) and o-(s) into the expression 
(f){sy / (f){s) = 7c{s)/a{s), the other part of the wave function is found as 

= (1 - qsy-^/'^s'^. (32) 

We write the upper spinor component in terms of the Jacobi polynomials 

<P+{s) = - s)^-^Pf v^' - 2s), (33) 

where -B„ is a normalization constant. The lower component of the spinor wave function 
can also be obtained by substituting Eq.(|H^ into Eq.(j2I). We should then solve the 
following equation 



a 



^ h ds 



[s . 



(34) 



where 7^ C . This is possible if a new variable is introduced as x = 1 — 2s . Now, the 
equation of the lower spinor component has been transformed into the following form 



a 



• ^ VoC{l-x) ^ jl-x) d 



e(l + x) 



dx 



(35) 



with 



0+(s) = C„(l - x)*v^(l + xf-'Pl^^^" ^\x), (36) 

where C„ is the normalization constant and its value is equal to i?„2'''^^~*^. If 
the following recursion relations and the differential formula satisfied by the Jacobi 
polynomials [TH] are included to the solution 

2 



;i + a;)P^ 



X 



X 



x) 



XI 



2n + jj, + Q + 1 
2 

2n + fj, + g + 1 



{n + g)Pi^' ^-^)(a:) + {n + l)Pitr'\x) 



x\ 



dx 



[X 



-n \ X + 



0- H 



n + fi)Pi^-'' ^')(x)-(n + l)P^^+i 

pi,, + g (n + /^)(^ + g) 



2n + jj + g 



pif, = ^ + + p(M. .+1), 



x] 



n + fi 



2n + n + g 



-pi'-r'\x) 



n-l 



x), 



(37) 



2n + iJ, + g + l " ' ' 2n + iJ, + g + l 

we obtained the lower spinor component in terms of the Jacobi polynomials as a function 
of s(r) 

;i - s)^s'^{{L + n/b)[{n + fi)Pj^^'^-^^ ^)(1 - 2s) 



")(1 - 2s)] +2{n + g)i2n + ^l + g + l- M)P^^'^'^ - 2s) 



-(1/26 +{n + 2h){2n + ^ + g + l)/2)P^2iVi, u)^^ _ 2^)} 
where L = CVq/^ + (1 + 7)/6 and M = ^ + iy^/b. 
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4. Conclusion 

We have solved the Dirac equation for the special form of the Woods-Saxon potential 
following a formahsm introduced by Alhaidari. The NU method is used to obtain a sys- 
tematical solution in the Dirac- Woods-axon problem. The energy spectrum of the bound 
states is analytically obtained and two-component spinor eigenfunctions are written in 
terms of the Jacobi polynomials. It is seen that the energy eigenvalues are a function 
of the parameter q and the solution space splits into two distinct subspaccs. We have 
seen that the non-relativistic limit of the Dirac equation can be obtained easily. We can 
also say that the exact results obtained for a special form of the Woods-Saxon potential 
give us some interesting applications in the various quantum mechanical studies and the 
relativistic nuclear scattering problems. 
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